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^ ' Abstract. The interest in decoherence-free, or noiseless subsystems (DFS/NSs) of 

^^, quantum systems is both of fundamental and practical interest. Understanding the 

invariance of a set of states under certain transformations is mutually associated with a 
^S] ■ better understanding of some fundamental aspects of quantum mechanics as well as the 

^ I practical utility of invariant subsystems. For example, DFS/NSs are potentially useful 

for protecting quantum information in quantum cryptography and quantum computing 
as well as enabling universal computation. Here we discuss transformations which are 
compatible with a DFS/NS that is composed of d-state systems which protect against 
[~^ ■ collective noise. They are compatible in the sense that they do not take the logical 

^D I (encoded) states outside of the DFS/NS during the transformation. Furthermore, it is 

shown that the Hamiltonian evolutions derived here can be used to perform universal 
quantum computation on a three qudit DFS/NS. Many of the methods used in our 
derivations are directly applicable to a large variety of DFS/NSs. More generally, we 
J^ ' may also state that these transformations arc compatible with collective motions. 
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1. Introduction 

In 1954, Dicke argued that an independent approximation for radiating physical systems 
was very often not reasonable [T] . Many physical systems are correlated, and transform 
in a similar, or even identical way so that an independent approximation should not be 
used. For some systems the motions can be considered collective. For example, this can 
be the case when a set of particles interacts with a common field. In proposed quantum 
computing devices, interactions with a common field can lead to collective noise which 
is an unwanted interaction. This collective noise, or collective decoherence, will cause 
quantum information to be lost to the environment [21 [3] . 

To avoid the detrimental effects of such noise, a theory of decoherence-free or 
noiseless subsystems (DFS/NSs) was formulated [H El El [TJ [U [9]. (For reviews see 
[TOl [TT].) When information is encoded into one particular type of DFS/NS it is 
protected against collective errors. More generally, under suitable conditions a DFS/NS 
can also protect against noises which are not of the collective type, but still correspond 
to an identifiable symmetry in the system-bath interaction. However, in this article 
we will focus on collective DFS/NSs. Once information is encoded into a DFS/NS, it 
must then be manipulated if the quantum information is to be used for the purposes of 
computation or simulation. Not just any physically available operation is acceptable. 
The manipulations, or gating operations, must be compatible with the DFS/NS ^ if the 
information is to remain protected. Here compatible means that the operations should 
not take the information outside of the DFS/NS. For, if they did, the information would 
be vulnerable to collective noises during the time it is not confined to the DFS/NS. 

The objective of this paper is to describe Hamiltonian evolutions which are 
compatible with collective operations on the system in the sense that they commute 
with the collective operations. The Hamiltonians we find here will provide a set of 
transformations which are compatible with a DFS/NS that protects information from 
collective noise and is composed of qudits. 

Qudits are quite interesting systems for several reasons. For example, two 
three-state systems, or qutrits, can be more entangled than two qubits [121 [131 [H]- 
There is also evidence suggesting that a collection of d-state systems can share a 
larger pairwise fraction of their entanglement capacity as the dimension increases |15j . 
They are beneficial for several information processing tasks including cryptography 
[111 m m [ini ED], computing [211 [221 [23], and games. They seem to be required 
for a version of the Byzantine agreement problem [23]. So for these reasons, among 
others, universality requirements for qudits [23 [26l (23, [2H1 [221 [HH] as well as several 
recent experiments and proposed experiments have provided methods for producing 
and manipulating qudits [311 1321 [331 [331 and an analysis of the collective behavior for 
an assembly of qutrits has been given [35] . 

Collective DFS/NSs have been studied quite thoroughly for a variety of reasons. 
One is that collective errors are physically observed in some systems. Another is 
that a collective DFS encoding can enable universal computation on a set of qubits 
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ElEZlESllMlliniSIllllllSlISSlllll- Collective DFS/NSs have also been observed 
to reduce noise in several experiments [l5l HHl HTl HHl US], including computation in 
a qubit DFS [501 [5T] . Even if collective errors are not present in an experiment, it is 
possible to induce such a symmetry using decoupling operations [52l [531 EH ESI ESI EI] ■ 
Theoretically collective noises are easier to treat in part because a basis for collective 
operations forms a representation of the algebra of the special unitary group [58] acting 
on the constituents. This enables a variety of group-theoretical methods to be employed 
in their treatment. We will see here that this is also the case for computing in a qudit 
DFS/NS. 

Specifically in this paper, Section [2] contains a review of DFS/NS theory. We then 
provide a description of the compatibility condition and show how it can be used to find 
a set of Hamiltonians which enable compatible computation on a DFS/NS in Section [3l 
Section H] contains a set of compatible Hamiltonians for a system of three qudits. We 
provide a complete set for qutrits in Section 14.11 and then generalize to qudits. We then 
show how to analytically obtain several unitary transformations on encoded states using 
these results in Section[5l Section[6]extends this analysis to a system of n-qudits {n > 3), 
thereby showing that these compatible Hamiltonians can be used to manipulate encoded 
quDit states {D ^ d) of arbitrary dimension. We conclude in Section [7] with a discussion 
of our results and their implication for universality using a three qudit DFS/NS. Two 
appendices have also been included which contain some detailed calculations. 

2. Decoherence-Free/Noiseless Subsystems 

Here a brief review of DFS/NSs is provided. Although the description is, by now, more 
or less standard, we will primarily follow the notation of Ref. [H]. For further details, 
see [7] and/or [8]. We note that one may describe a DFS/NS in terms of a collection 
of operators appearing in a semi-group master equation [6j. We will, however, use the 
Hamiltonian description here. 

2.1. Subsystem Structure 

Consider a Hamiltonian which includes a term describing the evolution of the system 
{Hs), a term describing the evolution of a bath or environment (Hb), and an interaction 
term describing the evolution of the system and environment together, (Hsb) 

H = Hs + Hb + Hsb. (1) 

Without loss of generality, we will write the interaction term as 

HsB = J2^'-^^-, (2) 

a 

where the Sa act only on the system and the Ba act only on the bath. Let us 
now consider the algebra A generated by the set of operators {Hs,Sa}- We will 
assume that this algebra is reducible and thus the operators can all be simultaneously 
block-diagonalized with a single unitary transformation. Such algebras are said to 
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be, mathematically, completely reducible, meaning they may be completely block 
diagonalized, i.e., written as a direct sum of irreducible components. It is this 
reducibility, which follows for any system-bath interaction symmetry that is preseved 
in the algebra, that allows for a non-trivial decomposition of the algebra into regions 
supporting the logical encoding of protected information. It is noteworthy that 
the "unitary trick" can be used to relate algebraic representation theory and group 
representation theory [58]. This is particularly useful for collective errors acting on 
systems of qudits since these operations form a representation of the algebra of such 
systems. 

This decomposition is described by the equation 

A=^lnj^Midj,C), (3) 

where the ?7,j-fold degenerate dj x dj complex matrices ^A{dJ,C) correspond to the 
irreducible components of A. We label these components by J, which collectively form 
the finite set JT". (It is important to note that this J actually stands for a set of 
quantum numbers when the constituents are qudits, with d > 3.) One may also define 
the commutant which is the set of elements that commute with all elements of the 
algebra A. The form of these is dictated by the irreducibility of the blocks in A, 

A' = {X:[X,A] = 0,\/AeA}. (4) 

To Eq. ([3]) there is a corresponding decomposition of the Hilbert space Hs = J2j ^"'^ ® 
C^'^ where the second factor corresponds to the part of the Hilbert space which is affected 
by noise (C'^-') and the first factor corresponds to that part which is not (C"-'). 

The unitary transformation which is used to to change between the physical and 
logical bases can be referred to as the DFS/NS transformation. In the logical basis, the 
elements of A exhibit a structure that allows quantum information to remain confined 
to the logical subspaces while the physical system interacts with its environment. The 
superpositions of physical states which form the logical, or encoded, states of the 
DFS/NS are created by this transformation in the following way. Let Vdfs be the 
aforementioned transformation that simultaneously block diagonalizes each element of 
the algebra Ai G A in an identical way, i.e., A[ = V^f^AiV^^. The physical states |\E'p) 
are then related to the logical states |\E'i) by 

|^L) = Kifs|^p). (5) 

In practice, Vdfs can be used to express states and operators in the physical bases in 
terms of states and operators in the logical basis as done explicitly, for example, in |62j . 
We may now define a decoherence-free or noiseless subsystem in the following way. 
Suppose that we represent a basis of eigenstates corresponding to a particular J with 
the set {|A) (8> |/i)}, where A = 1, ...,nj and yU = 1, ..., dj. Then if 

dj 

Aa|A)®|/i) = 5^M^^,,,|A)®|/i') (6) 
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for all Aa, A, and /x, there exists an irreducible decomposition as given by Eq. ([3]). The 
invariance of the degeneracy labels |A) in this last expression reflects the ability to 
reliably store quantum information in certain regions of the system Hilbert space when 
the algebra A can be decomposed in this way. The information is stored in blocks with 
the same J but different A. Each A specifies a particular DFS/NS basis state. These 
logical states can therefore be expanded in terms of those states associated with a given 
A. Although the initial encoding of a particular logical state may change, it will remain 
confined to its initial subspace. 

A decoherence-free subspace is one for whcih the matrices M are numbers (1x1 
matrices) which act on a one- dimensional representation, i.e., a singlet state. 

2.2. DPS and NS examples 

In this section we review a few examples of DFS/NSs which will aid in our discussion 
of qudit systems. The examples will be useful for comparing and contrasting certain 
properties of qubit and qudit systems. 

2.2.1. Four-qubit DFS As stated above a decoherence-free subspace is comprised of 
singlet states. Four qubits can be used to construct a DFS/NS qubit which is represented 
by two singlet states, one singlet state for the logical zero and one for the logical one 
[6]. The logical states considered here protect the information from errors which act 
the same on each of the four physical qubits constituting the system. These states are 
given explicitly by 

|0l) = (10101) + 11010) - lOllO) - |1001))/2, (7) 

and 

\Il) = (2 |0011) + 2 |1100) - |0110) 

- |1001) - |0101) - |1010))/v^, (8) 

where |0) and |1) represent two orthogonal basis states for a spin-1/2 particle. 

It turns out that the Heisenberg exchange interaction is universal for a set of qubits 
constructed in this way and logical gates consisting of only these interactions have been 
provided [36l El ETJ [59] . (In this paper we neglect the corrections which must be made 
to the logical gates in order to account for three- and four-body interaction terms. 
For a discussion of these effects see, for example, Ref. [60].) The Heisenberg exchange 
interaction between pairs of physical qubits can be expressed as: 

E,, = -{I + a,-a,), (9) 

where a = {crx,cry,az) is the vector of Pauli matrices. As written Eij is the exchange 
operation between qubits i and j, i.e., Eij\(j)) i\ip) j = \ip)i\(j))j. The logical "X" operation 
is given by 

X = ^{E2s-E,,). (10) 
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The logical "Z" operation is given by 

Z = -Eu (11) 

and Y can be obtained from these two by commutation or the finite transformations 
can be written in terms of Euler angles. 

2.2.2. Three-qubit DPS The smallest number of qubits that enables the DFS/NS 
encoding of a logical qubit is three. In this case the logical zero and logical one are 
each represented by a doublet. In the {|0) , |1)} basis, these states can be written as 

|0i) = ao(lOlO) - |100))/V2 + /5o(|011) - |101))/V2, (12) 

and 

|1l) = ai(2 |001) - |010) - |100))/v^ 

+ /5i(-2|110) + |011) + |101))/V6, (13) 

with |ajp + |/5jp = 1. As in the previous example, these logical states protect 
the quantum information from collective errors. Although the Heisenberg exchange 
interaction is universal for both the three-qubit and four-qubit DFS/NSs, there are 
some practical differences between the respective logical operations for these systems 
with regard to error prevention [611 [62]. We will not discuss this here, but will note the 
form of the logical operations. The logical "X" operation is given by [H] 

X = i=(E23-i?13), (14) 

while the logical "Z" operation is given by [63] 

Z=i(Ei3 + ^23-2^12). (15) 

Again, Y can be obtained from these two by commutation or Euler angles. 

It is important to note the similarities and differences here. However, for our 
purposes, the most important point is that the Heisenberg exchange interaction can be 
used to construct a universal set of operations for both the three- and four-qubit DFS. 

2.2.3. Three-qutrit NS Here we provide some details concerning the logical qubit 
encoding over a subspace of three physical qutrits. The purpose of this particular 
example is to provide enough structure to enable a smooth transition into our discussion 
of qudit systems. 

The three qutrit DFS/NS is quite analogous to the three qubit DFS/NS discussed 
above in many respects although the dimensions of the subspaces are different. For the 
three qubit NS, a tensor product of three qubits can be decomposed into two two-state 
subsystems and a four-state system. This may be written as2®2®2 = 2©2©4. A 
tensor product of three qutrits can be decomposed into a singlet, a decuplet, and two 
octets: 

3®3®3 = 8©8©1©10, (16) 
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where the 8 identifies an eight-dimensional representation (octet), the 1 identifies a 
one-dimensional representation (singlet), etc. A logical qubit can be represented by 
two degenerate eight-state subsystems, the two octets [58]. We may also note that it 
is possible to experimentally produce such states encoded in a set of polarized photons 
[32]. 

Exphcit forms for the states within the two octet subspaces of three physical qutrits 
were provided in [58]. The three states within each individual qutrit will be denoted 
|0), |1), or 1 2). In terms of the quantum numbers p, q, t, t^, and y as given in [58], these 
are 



(17) 



|0) = 10,1,-1/2,-1/2,-1/3), 
|1) = 10,1,-1/2,1/2,-1/3), 
|2) = 10,1,0,0,2/3). 

The physical, or computational basis (|000), |001), ... , |222)) transforms into the DPS 
basis via the operation Vdfs which is the matrix composed of the Wigner-Clebsh-Gordan 
coefficients. In this noiseless basis the logical zero state \0l) is formed from the octet 
which corresponds to the states [58] 



8,0 _ 
1 " 


= (|200) - 


|020))/v^. 


r- 


= ( 100) - 


|010))/v^. 


r- 


= (011)- 


|101))/v^. 


r = 


= (211)- 


|121))/v^. 


r- 


= ( 122) - 


|212))/v^. 


r- 


= ( 022) - 


|202))/V2, 


8,0^ 


= (-1021) 


- 120) + 1 



?/'8'° = (2|012) + |021)-2|102) 

- |120)- |201) + |210))/v^, 

where the first superscript on each ip denotes the dimension of the representation, the 
second is a degeneracy label and the subscript labels the state within the representation. 
A second octet forms the logical one state \Il), 



r = 


;-2|002) + |020)+ 200))/y6 


'/ = 


;-2|ooi) + |oio)+ ioo))/Vq 


r = 


-2|110) + |011)+ 101))/V6 


r = 


-2|112) + |121)+ 211))/V6 


f = 


-2|221) + |122)+ 212))/V6 


r = 


;-2|220) + |022)+ 202))/V6 


v = 


;-2|012) + |021)-2|102) 



(19) 



^ 



+ |120) + |201) + |210))/V12, 
(|021)- |120) + |201) - |210))/2. 
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In terms of these octets the logical zero state is given by an arbitrary superposition 
of the eight ^jJj' states, \0l) = Tlij^ji^j' ^^"^ likewise for 11^;,) = ^jPjip/ ■ According 
to the theory of noiseless subsystems, states within the sets {ipi' } ({V'i' }) will mix 
together in the presence of collective noise but not with the states in {ijj^' } {{ip^' }) 
or with the one- or ten-dimensional representations. Furthermore, the mixing will be 
identical for both octets. This means that when only collective errors are present, the 
information encoded in lipi) = a \0l) + b \1l) will be protected. It is also important to 
note that the initialization is arbitrary in the sense that an arbitrary combination can 
be taken [65]. In practical situations this can be very beneficial. (See for example |32j.) 

2.2.4- Three-qudit DPS As discussed in [5H] the tensor product of three equivalent 
irreducible representations of SU((i) gives rise to the smallest number of qudits for which 
a NS, representing a qubit in terms of qudits, exists. This can be seen in the tensor 
product of three d-state systems using the tableau of a representation. Taking the tensor 
product of three such systems produces two degenerate tableau. These correspond to 
two degenerate irreducible representations which can be used to store protected quantum 
information in the form of a noiseless qubit. 

As we will see later, this decomposition is valid for all three-qudit systems and 
enables us to calculate the compatible transformations governing the noiseless evolution 
of such systems. 

3. Computing in a DFS/NS 

When information is encoded in a DFS/NS, it will be protected during the unitary 
transformations which are used to manipulate it only if those transformations do not 
couple states inside the subspace with states outside of the subspace. In other words, 
we must restrict the transformations to those which preserve the subspace structure. 
When this can be accomplished, we say that such logical operations are compatible with 
the DFS/NS. In order to find a DFS-compatible set of gates we will find it useful to 
first identify the stabilizer of a DFS. 

3.1. The Stabilizer of a DFS/NS 

Let |\E') G C be a state in the DFS, where C denotes the code space. Let us define a 
stabilizer, iS, as in [8]: 

S = {S: S\^) = |^),V|^) eC}. 

While the elements of this set do in fact leave all of the code words unchanged, we can 
relax this requirement slightly when considering the DFS states since in this case the 
information being encoded is stored in the states labeled by A. This property of the 
DFS/NS states allows us to define a modified version of the stabilizer above as 

dj 

S' = {S' : S'{\X) ® l/i)) = \X)^Y1 ^/^A^' I-"') '^ 1"^) e ^}- 
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We can parametrize elements S' of the modified stabilizer which are relevant to the 
system using Eq. ([6]) and a set of arbitrary complex numbers {fa}: 



D{vi,V2,...) =exp 



"^VaAc 



(20) 



This leaves the labels A unchanged, as seen by expanding the exponential and acting 
term by term on states of the form |A) ® |/i). 

3.2. Compatibility Conditions for DFS/NS Evolution 

Using the modified stabilizer, we can state that U is compatible with C, if V|\E') G 
C, f/(|A) ® l/i)) = Ea'^aa'IA') ® E;,'^/./.'!/^') and therefore S'U{\X) ® |^)) = 
S'iEx' -hw |A') ® J:,' K,,. I/x')) = (Ev -hw |A') ® E^' K'^^' I/"')) VS' G S' . This implies 
that U-^S'U{\\) ® 1^)) = |A) ® Y.^' Kt^' \l^') so that 

U~^S'U e S'. (21) 

The condition Eq. fl2Tl) may now be re-expressed as 

UD{vuV2,...)U^ = S', (22) 

for some S" G iS', or 



C/exp 



^f^Ao 



f/t 



exp 



Y^VaUA^U^ 



S'. 



Taking the natural logarithm of this equation gives 
Y^VaUA^U^ = ln{S'). 



(23) 



(24) 



Now let us define U = exp [—iHt]. Then, after taking the derivative of both sides of 
this last expression with respect to time, we find a very simple and sufficient, but not 
necessary condition for the compatibility of Hamiltonians 

[H,A^] = 0,^A^eS'. (25) 

Then U = exp{—iHt) is compatible with the DFS/NS. Clearly if the Hamiltonian 
which generates the unitary U commutes with every element of the algebra, then U 
also commutes with every element of the algebra. This provides a subset of the unitary 
transformations U satisfying Eq. (I2T!) . 



4. Explicit forms for Compatible Hamiltonians 

It was shown in the previous section how a parametrization of the stabilizer in terms of 
time independent coefficients f „ can be used to identify a set of unitary transformations 
that are compatible with the DFS/NS states. This section provides a set of those 
Hamiltonians satisfying Eq. ( l25l) for systems composed of three qudits. In the following 
analysis we neglect the internal Hamiltonian of the system and assume that the only 
interaction terms present are those which can be considered as collective. The set 
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presented here is valid for all d > 3 and it is unique for the case oi d = 3. The algorithm 
for determining this set could be used for any DFS/NS, not only those which protect 
against collective errors. 

4-1- Explicit forms for Qutrits 

Using MATHEMATICA, we have determined the complete set of Hamiltonians which are 
compatible with an encoded two-state noiseless subsystem of three qutrits (compatible in 
the sense that they satisfy Eq. fl25l) ). To do this, we expanded the relevant quantities in a 
complete set of traceless, Hermitian matrices {Aj}. For these calculations, we employed 
the Gell-Mann matrices, but in principle any basis can be used. Several properties and 



conventions of these matrices are given in Appendix A 



Collective errors will be denoted Sa-, of which there are eight. For example. 
Si = ^j A^ , where the superscript labels a particular qutrit and the subscript denotes 
the type of error labeled 1 through 8. Let Hijk = \i \j ® A^, where the Xj are 
the Gell-Mann matrices with Aq = 1. An arbitrary Hamiltonian can be expanded as 
^ — Tliijk (^ijkfJ'ijk with H traceless if not all i,j, k are simultaneously zero and the aijk 
are arbitrary real constants. The algorithm we used to find the H which commutes 
with every collective error proceeds as follows. (We emphasize that this algorithm is 
quite general and could be used with the elements of any stabilizer to find compatible 
Hamiltonians.) 

(1) Expand H in terms of the complete set of Hermitian matrices ^ijk described above: 

(2) Determine the commutator of the general Hamiltonian H and a generic collective 
error Ylidi^i^ where the Qi are arbitrary coefficients. In other words, calculate 

(3) Find the projection of [H, ^^ QiSj] onto a component ^ijk by taking the trace 
of the basis element ^ijk with the result of (2). In other words, calculate 
tr ([if, Y.i 9iSi] fiijk) for each fiijk. 

(4) Set all of the projections equal to zero and then solve the system of linear equations 
for the expansion coefficients a^jk which satisfy these relations, thereby determining 
the H which will commute with Ylidi^i- 

The results are summarized in the following 

8 

hi = ^/ion, (26) 

i=l 
8 

h2 = ^/ijoi, (27) 

i=l 
8 

^3 = J^ IJ-iio, (28) 
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hi = 2_^ fijkfJ^ijk, (29) 

ijky^O 

^5 = 2_^ dijkfiijk, (30) 

where the fijk are the structure constants and the dijk are components of the totally 



symmetric d-tensor. (See Appendix A ) These are the only Hamiltonians, along with 
combinations of these, which commute with all of the collective errors that may act on 
a set of three qutrits. The first three are generalizations of the Heisenberg exchange 
operations which are known to be universal for the three- and four-qubit DFS/NSs [8]. 
The last two are, perhaps, not immediately obvious candidates for DFS compatible 
Hamiltonians. However, we will show analytically that these two, along with the first 
three, are compatible with all three qudit NSs. 

4-2. Hamiltonians Compatible with a 3- qudit DFS 

We will now show that a set of Hamiltonians having the form of Eqs. (I26l) - (l30l) are 
compatible with a three-qudit DFS/NS and can be obtained analytically. This follows 
from the fact that Eq. (1A.2I) . the Jacobi identity Eq. (1A.3I) . and Jacobi-like identity 
Eq. (jX4]) hold for all SU(ci), d > 3 [66]. 

Let us consider a Hamiltonian constructed from the fiijk with the constituent basis 
elements Aj belonging to a matrix representation of the Lie algebra of SU((i) with d 
arbitrary: 

H = 2_^ CLijkfXijk- (31) 

ijk 

The three qutrit case above is a special case for SU((i) when d = 3. If we consider 
collective errors on the DFS/NS, 

Sj = /ijOO + f^OjO + /^OOj) 

our condition, Eq. fl25l) reads 

[S,,H] = 0, 

which implies that we require 

film.O'jkl + fjlrndkil + fklmdijl = 0. (32) 

Comparing this equation with Eqs. ( JA.31) and f lA.4l) . we can immediately see that there 
are two sets of numbers aijk which satisfy Eq. (132|1 . fijk and dijk- In other words, this 
equation is satisfied for aijk = fijk and for aijk = dijk- 

To obtain Hamiltonians of the form hi,h2,h3, we let one of the factors be the 
identity in fiijk. Let us consider one particular case, the case when the third factor is 
the identity 

TV 
H = 2_^ O'ijOfJ'ijO 
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and calculate [H, Sj] . The result is 

O-mflkj + (J'ljoflki = 0. 

If we multiply by dmki and sum over k and i, we get 

(Imflkjdmki = 0. 

Clearly, if auo = Su this equation is satisfied and we obtain the form Eq. ( l28l) . Similarly 
for Eq. (!26l) and Eq. (l27j) . We emphasize that the analytic proof is valid for all three- 
qudit DFS/NSs, not just qutrit states. 

Also, we have not proved the converse; i.e., we have not shown that these are the 
only Hamiltonians which commute with collective errors on three qudits. The algorithm 
in the previous section has been implemented for qutrits and, in that case, these are 
the only Hamiltonians satisfying the stated conditions for DFS/NS compatibility. We 
suspect that this is also the case for qudits. 

4 ■ 3. Generalizations 

It is important to note that several of these results are quite general. First, we reiterate 
that the algorithm used to find the Hamiltonians for the collective three qutrit DFS/NS 
can be used for any DFS/NS given a basis for the stabilizer elements of the code space. 
Second, as shown later, due to the permutation symmetry of qudits undergoing collective 
errors, the Hamiltonians we have derived are compatible with any qudit DFS/NS 
undergoing collective errors. Perhaps even more generally stated, the d-state system 
Hamiltonians with the same form as Eqs. (!26|) - (l30|) . including the generalized exchange 
interaction, commute with any collective transformation on a set of qudits. 

Furthermore, as shown in the next section, the exchange Hamiltonians can be 
analytically exponentiated to provide generalized SWAP operations. The SWAP 
generates the permutation group on the set of qudits and this group commutes with the 
group of collective unitary transformations. A description of this is given in Ref. [30] 
where it is shown that one can use this relation to produce efficient qudit circuits. 

5. Universal Computation 

We may now ask the question, are these Hamiltonians sufficient to perform arbitrary 
unitary transformations on NS qubits? In other words, can we perform universal 
computation on this NS using the compatible Hamiltonians provided here? 

In order to answer these questions, we will first note some important properties of 
the algebra of the Hamiltonians we have found to be DFS-compatible. Then we attempt 
to find the corresponding unitary transformations. 
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5.1. Commutation Relations 

Now let N = d'^ - 1 be the number of matrices in a basis for the Lie algebra of SU{d). 
We will define the (i-dimensional analogs of Eqs. ( l26l) - (130|) by 

N 

ei(c?) = "^l^oii, (33) 

i=l 

N 

62 (C?) = 5]^/iiOi, (34) 

TV 

^3(c^) = X^/^no, (35) 

^(^) = 5Z fijkl^ijk, (36) 

^(^) = Xl ^iJkfJ'ijk- (37) 

where the /ijjA; now represent the tensor product of three d x d basis matrices. It can 
be shown that the commutation of any two different ei{d) Hamiltonians yields the F{d) 
Hamiltonian, 

[ei{d),e2{d)] = -2iF{d), (38) 

[ei{d),e3id)] = +2iF{d), (39) 

[e2(rf),e3(rf)] = -2iF{d). (40) 

One can also show that the commutation of an ei{d) Hamiltonian with the F{d) 
Hamiltonian gives a combination of the remaining two e^'s. 

[ei(rf),F(d)] =4^(e2(ci)-e3(rf)), (41) 

[e2{d),F{d)] =At{e^{d)-eM). (42) 

[eM,F{d)] =At{e,{d)-eM)- (43) 

Furthermore, we have found that the D{d) Hamiltonian commutes with the three 

Cj's as well as F{d), 

[aid), D{d)]= 0,1 = 1,2,3 (44) 

and 

[F{d),D{d)] = 0. (45) 

This last relation, Eq. (H5|) . can be obtained by an explicit expansion of the two ordered 
products using Eq. (]A.2p . The expansion may be reduced to terms involving products 



of the form fijmfkim which can be expanded using Eq. ( JA.91) . What remains can be 
reduced further with the help of Eqs. flA.12p and flA.13P giving the stated result. 

With these results, we may now show that a sub-algebra isomorphic to the Lie 
algebra of SU(2) is generated by a combination of these Hamiltonians. First, note that 
from Eqs. (14T|) and ( l42l) . we may show that 

[(ei-e2),F]=4z(ei + e2-2e3). (46) 
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From Eqs. ([MD-din]) 

[(ei - 62), (ei + 62 - 263)] = -l2tF. (47) 

Therefore, the three matrices (ci — e2)/2^/3, (ei + 62 — 2es)/6, and F/2\^ form a 
representation of the Lie algebra of SU(2). We will now use this result in the construction 
of the logical analogues of the Pauli matrices acting on the encoded qubit states. 

5. 2. Unitary Transformations and Logical Operations for Qudits 

Before obtaining the logical gating operations, it is interesting to note that analytic 
expressions for the exponential of each one of the three Cj Hamiltonians in the physical 
basis may be obtained for qudits. There are two fortuitous properties of these matrices 
which enable us to provide such an analytic expression: 1) the sum of the off-diagonal 
matrices commutes with the sum of the diagonal matrices and 2) when squared, the sum 
of the off-diagonal matrices is diagonal and can easily be summed. This can be shown 
to be true by direct computation using the set of Gell-Mann matrices in the case of 



qutrits. [Appendix B| provides proof that it is also true for qudits. These two properties 
enable us to sum the series resulting from the exponential of the Hamiltonians ei, 62, 63. 
An explicit form for the unitary evolution corresponding to these Hamiltonians is also 



given in [Appendix B 



Now let us define the logical "X" operator, X, which acts on the DFS through the 
relation 

X = ^(61-62). (48) 

We note that the overall sign of the states spanning logical one Eq. (TT9!) are chosen so 
that the form of X above resembles the expressions appearing in Eq. flT^ and Eq. 01U|) 
for the three- and four- qubit DFSs, respectively. It should also be mentioned that 
we are describing the logical X operation in terms of Hamiltonians rather than unitary 
transformations. For comparison, notice that Eqs. (IT^ and (ITU]) could also be written 

as ^(<?2 ■ 0^3-^1 •(T3). 

The exponentiation of X leads to a time evolution given by 

?7x = ll + iXsin(t)-X2(l-cos(t)). (49) 

This can be obtained in two different ways. One way is to use the prescription for the 
exponential of the Cj as described above. Another one is to calculate the exponential 
from the results of Appendix B.3 and then transform back to the physical basis using 



the fact that for a similarity transformation V , and a Hamiltonian if, 

V Q^^{-iHt)V-^ = exp{-iVHV-H). (50) 

This enables one to transform between the physical, or computational basis states and 
the logical basis states as well as between the different sets of operators - the physical 
unitary transformations one would implement in experiments and the logical ones. 
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Similarly, the logical "Z" operator, Z, can be expanded in terms of the e^'s by 

Z = 1(61+62 -263). (51) 

o 
The form of the unitary transformation resulting from the exponentiation of Z is similar 
to that of Ux and is given by 

[/z = l-«Zsin(t)-Z2(l-cos(t)). (52) 

Once X and Z are found the logical Y can be produced through commutation. 
This shows, with our particular choice of scaling, that the three matrices proportional 
to (ei + 62 — 263), (e2 — Ci) and F form a representation of the Lie algebra of SU(2) and 
act as the Pauli matrices on the DPS qubit made from three qutrits. 

A rotation in SU{2) about an arbitrary axis can be obtained by three successive 
rotations. In particular, a rotation about the logical Y may be performed using the 
following decomposition of the logical SU{2) group, using real (Euler) angles a,/3, and 

7: 

U{a,P,'y) = exp[—iZa]exp[—iX.p]exp[—iZ'j]. (53) 

Indeed this parametrizes all of the group SU(2). Therefore, we have shown that for an 
encoded qubit comprised of three qudits, these DFS/NS compatible operations alone 
can perform any rotation over the logical subsystems. 

5. 3. Unitary Transformations and Logical Operations for Qutrits 

In the logical basis, f/x S'Cts as the identity on the singlet and decuplet states (given in 
Eq. ( TT6l) . and explicitly in [5H]) throughout the entire evolution. The action of X on the 
states forming the logical zero Eqs. (fT8|) is such that it swaps the states for their logical 
one counterparts Eqs. ( iT9l) and vice versa, 

X^r = ^?\ (54) 

X^r = ^r • (55) 

In other words, it acts as a Pauli X gate on the logical states. These relations can be 
used along with Eq. (H9l) to calculate the action of ?7x o^ logical zero basis states, 

Ul^ii^T) = [^ + ^Xsin(t) - %\l - cos(t))]^J'° 
= i>j' +i sm{t)'iljj' - (1 - cos{t))ipj' 
= cos(t)ipj' + ism{t)ipj' . (56) 

It therefore follows that 

t/x |0l) = cos(t) \0l) + i sin(t) \U) , (57) 

and similarly, 

?7x I 1l) = cos(t) I U) + 2 sin(t) \0l) . (58) 

Now note that the logical "Z" operator, Z, acts as the identity on states in octet 1 while 
changing the overall sign of states in octet 2 and this can be used to obtain Eq. (ISTIl 
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directly for qutrits using the explicit expressions. The unitary may then also be obtained 
directly from Z. The transformation of logical zero is given by 

Uz |0l) =U2J2 '^J^f'T = [1 - ^Z sin(t) - 2^(1 - cos(t))] ^ a,V|'° 



y ajijjj' — i sin(t) \^ ctji^j' — (1 — cos(t)) y^ otjil^j' 



(59) 



= |0L)exp(-2t), 
while logical one transforms unitarily by 

f/z|U) = |lL)exp(+^t). (60) 

Again, the decuplet states are left unchanged by the action of f/z. This implies, along 
with the invariance of the singlet state, that this gate set is canonical in the sense 
described in [62] which is important for applications of decoupling pulses to eliminate 
leakage and protect the information [611 [62] . 



5.4- SWAP Operation 

The exchange, or SWAP operation between qudits p and q can be achieved by allowing 
the appropriate DFS/NS compatible Hamiltonian em{,d), {p ^ m ^ q) to act between 
the two for a specific amount of time. To show this, let us write 

^ A, ® A, = ^Ai® Ai + ^Aj ® Aj, (61) 

s=l i j 

where the Aj (Aj) represent the diagonal (off-diagonal) components of the traceless, 
Hermitian basis {A^} normalized such that Tr(AsAs/) = 26ss'- {In this section only, we 
use i for indices on diagonal elements of the algebra and j for elements of the algebra 
which have no nonzero diagonal elements.) Since ^^ Aj ® Aj commutes with ^ . Aj ® \j 
(see Appendix B.l ) we may express the exponential of X]s=i ^s ® A^ as 



exp 



-it 



it ^ A^ 



A. 



s=l 



exp 



-it 



^ Ai (g) Ai 



exp 



it 2_^ Aj ® Aj 



-it 



(62) 



The off-diagonal members of this basis may be written either as \k){l\ + \l){k\ or 
i \k){l\ — i \l){k\ for k ^ I (for notational simplicity we use the two-qudit computational 
basis ({|11) , 1 12) , . . . , \dd)}), so that the off-diagonal contribution appearing in this 
last equation may be written as 



exp 



—it 2_] Aj ® Aj 



exp 



-it^Mk. 

k<l 



(63) 



where 



Mu,i = {\k){l\ + \l){k\) ® (|A;)(/| + \l){k\) 

+ ii\k){l\-t\l){k\)^{t\k){l\-t\l){k\) 
= 2\k){l\(g)\l){k\+2\l){k\(g)\k){l\. 



(64) 
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For SU{d) the number of M^^i appearing in the summation of Eq. fl63p is (rf^ — d)/2. 
Also, each distinct M^^i can be seen to commute with the others since {\k){l\ (8> |0(^l + 
\l){k\ ® \k){l\) X {\p){q\ (g) \q){p\ + \q){p\ (g) \p){q\) = when p ^ I jL q a.nd p j^ k j^ q. This 
allows us to write 



exp 



-it 



J2 ^j ® ^j 



exp [~itMi^2] ... exp [-itM^-i^] . (65) 



When squared, each M^^i is diagonal and given by M|; = 4 \k){k\ ® |/)(/| + 4 |/)(/| (g 
\k){k\ , which implies that M|; = 4Mfc^/, etc. Now, if we choose to define 

Qk,i = ^ = \k){l\ ® |/XA;| + |/)(A;| ® |A;X/| , l^k (66) 

and 

Rk,i = ^ = \k){k\ m\ + m\ \k){k\ , l^k (67) 

we may express the exponential of M^^i as 

Uk,iit) = exp[-^tMfc,;] = 1- tQk,i sin(2t) + /?fc/cos(2t) - 1). (68) 

If we now let t = tt/A we find that 

Uk,i{t = 7r/A) = l-tQk,i-Rk,i. (69) 

Using the definition provided in Eq. fl66l) we see that the only two-qudit computational 
basis states |a/3) which survive the action of a particular Q^^i are \lk) and \kl), i.e., 

Qk,l \aP) = Sa,l5f3,k \kl) + (5a,fc5/3,i \lk) . (70) 

Also, using Eq. (!67|) . we find that 

(1 - Rk,i) \a(3) = \aP) - Sa,iSi3,k \aP) - Sa,kSp,i \aP) . (71) 

Therefore, at t = 7r/4 the two-qudit computational basis states |a/5) evolve to 

t/fc,Kvr/4) \af3) = \af3) - 6^,i6p,k{t \kl) + \a(3)) - <5,,fc<5^,/(* 1^^) + W(^))- (72) 

In other words, Uk^i{TT/4:) acts as the identity on all product states \a(3) except \kl) 
and \lk). On these states ?7fci(7r/4) exchanges \kl) {\lk)) for |/A;) {\kl)) along with a 
phase shift of exp(— 'J7r/2). The action of Gxp[—itJ2s=i -^^ ® K] at t = 7r/4 on a given 
product state |a/3), a, P = 1,2, . . . d can now be calculated using Eqs. fl^ and 1^^. 

<f-i 



exp 



/4) J] A, ® A, 



s=l 



. exp(-2(7r/4)X;iAi(g)Ai) |aa), if a = /3, 
"^/ ^ _, exp(-«(7r/4) E. A. ® K) |/?«) , if a ^ /?. ^ ^ 



In order to determine the diagonal elements of exp(— it^^ Aj ® Aj) we must first find 
the coefficients ^^.n which satisfy the relation 

d 

^ Aj (g) Aj = ^ i^^n \m){m\ ® \n){n\ . (74) 



m,n=l 
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In the following analysis we will take a Gell-Mann basis for SU{d), in which there 
are a total of c? — 1 diagonal elements. All of these diagonal matrices can be constructed 
using Eq. ( JB.81) . For SU{3) the two diagonal elements are given by A3 = |1)(1| — |2)(2| 
and As = (|lXl| + |2)(2| -2 |3)(3|)/V3. Two of the three diagonal elements of SU{A) can 
be formed by respectively placing A3 and As in the upper left block of a 4 x 4 matrix 
with zero entries on the fourth row and column. The remaining diagonal element is 
then given by Eq. (^M) with rf = 4, i.e., A15 = (|1)(1| + |2)(2| + |3)(3| - 3 |4)(4|)/v^. 



This procedure can be used to construct the set of diagonal matrices corresponding to 
SU{d) for all d > 3; one simply places the diagonal elements associated with SU{d — 1) 
into the upper left block oi a, d x d matrix with zeros along the d*^ row and column. 
Eq. (IB. 81) can be used to obtain the remaining element. 

As we proceed to determine the expansion coefficients ^rn,n appearing in Eq. ( !74l) 
let us now write the sum of off-diagonal matrices explicitly as 

5^A,®A, = (|1)(1|-|2)(2|)®(|1X1|-|2X2|)+ ... 

i 

+ ^(^^(|iXi| + ■■■ +\d- i){d -i\-{d-i) \d){d\) 

® (|1X1| + ... +\d- l){d - 1| - (rf - 1) I^X^I). (75) 

If we consider the case when m = n we find that the coefficient C,m,m attached to the 
matrix element \m){m\ ® \m){m\ obeys the relation 

{ \i , a m = d, 

2(m—l) I n sr^d, 1 -r / i \'^) 

— + 2 > ,1 , ,, , it m jt d. 

Let us now evaluate the series which appears in the expression for C,m,m when m ^ d. 

A 1 _ 1 1 1 1 

^£r+i -"(-" - 1) ~ '^('^ + 1) ^ (m + l)(m + 2) + ••• +(rf-2)(d-l) + (rf-l)d 



111 111 



m m + 1 771+1 d — 1 d— 1 d 

1 1 d — m 



m d rad 
Therefore, when m ^ d we find that the coefficient C,m,m can be expressed as 



(77) 



so that 



2(m-l) ^ ^ y 1 _ 2(m-l) ^ 2{d - m) _ 2{d - 1) ^ 

m -^^ u(u — 1) m md d 

lj.=m+l '^^'^ ' 

em,m = ^^^^, form = l,2, ...,c?. (79) 

d 

Now suppose that m ^ n. Examination of the expansion given by Eq. fl75l) reveals 

the symmetry of all terms \m){m\ ® |'^)(^| under the interchange of indices m and n. 

Since the relation ^rn,n = in,m then follows, we may assume that n > m in the proceeding 
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discussion without loss of generality. Given this assumption, we find that the coefficient 
Cm,n appearing in Eq. fl7i|) satisfies the equation 



^m,r, 



_2 



if n = d, 
if n ^ d. 



(80) 



However, this relation for ^rn,n when n ^ d may be simplified in view of Eq. (I77|l . 

d 



- + 2E 



tJ.=n+ 



1 _ 2 2{d-n) _ 2 
fi{fi — 1) n nd d 



And so we find that all of the coefficients ^m.n {fn y^ n) are equal and given by 

2 



<,m,r, 



d' 



for m,n = 1,2, . . . ,d and m ^ n. 



12) 



The exponential of ^^ Aj ® Aj = ^^ „^^ ^m,n |m)(m| ® |n)(n| can now be evaluated. 



exp 



-it ^ Aj (g) Aj 



2J exp(-it^m,„) |m)(m| (g) \n){n\ 

m,n=l 

y ^exp(— (2it((j — l))/(i) |m)(r7i| (8> |m)(m| 



+ Y^ exp(2it/(i) |m)(m| ® |n)(n| . (83) 

Therefore, the unitary transformation which corresponds to the exponentiation of 
^5=1 As ® As can be expressed using Eqs. (1^ . (1^ . (1^ . and (15^ . In particular, the 
two-qudit computational basis states |a/5), a, (3 = 1,2, . . . d evolve at t = 7r/4 to 



exp 



-i(7r/4) ^ As ® As 



\a(3) 



exp(—(7Ti(d—l))/2d)\aa), if a = /3, 



-iexp{7Ti/2d) \i3a) 



a ay^(3. 



U) 



However, since exp(— (7rz((i — l))/2d) = exp{—'Ki/2)exp{ni/2d) = —iexp{'7ii/2d), we 
have 



exp 



-2 vr 



/4) ^ As ® As 



s=l 



\aP) = —iexp{Tci/2d) \Pa) for a, 13 = 1,2, . . . ,d. (85) 



^d2_l 



Although the exponential of ^s=i -^'^ ® -^s induces a phase shift on all product states 
\aj3), the shift is the same for all initial states as they evolve to t = 7r/4. This unitary 
transformation thus acts as the SWAP operation at t = 7r/4 between any two d-state 
systems. 



6. Hamiltonians Compatible With an n-Qudit DFS 

In this section we show that the Hamiltonians given by Eqs. fl33|) - fl37|) remain compatible 
with a collective DFS/NS encoding of n physical qudits into a set of logical quDits, 
d ^ D. Since, in principle, a logical quDit of dimension D can be formed by such 
an encoding, the Hamiltonians considered here have the ability to generate compatible 
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transformations on an encoded quDit state of arbitrary dimension. To see this, let us 
first consider the exchange operations ei{d) which were previously shown to commute 
with the collective errors, i.e., 

[ei(rf), Sj] = [1 O Aj (g) Aj , Aj (g) 1 (g) 1 + 1 O Aj (g) 1 + 11 O 1 (g) Xj] 
= 1 (g (AjAj - XjXi) (g Aj + 11 (g Aj (g {XiXj - XjX,i) 
= 2ifijk{l ® Afc ® Ai - 11 ® Afe ® Ai) 
= 0. (86) 

(Recall that the structure constants fijk are totally antisymmetric and that 
the generators of SU{d) obey the commutation relations [Aj, Xj] = 2ifijkXk- 
Again, summation over repeated indices is implied.) Notice how the term 
[Aj (g Aj (g 1 , 11 (g 1 (g Xj] in the commutator of [ei{d), Sj] can be immediately neglected 
while the nontrivial contributions are proportional to [Aj, Xj] (g Aj + Aj (g [Aj, Aj]. The 
commutators of e2{d) and 63(6?) with the collective errors also share these properties, 

[e2{d), Sj] = [Aj (g 11 (g Ai , Aj (g 11 (g 11 + 11 (g Aj (g 11 + 11 (g 11 (g Xj] 
= [Xi, Aj-] ® 1 ® Ai + Ai ® 11 ® [Xi, Xj] 
= 2ifijk{Xk ®t®Xi-Xk®l®Xi) 
= 0, (87) 

and similarly for [e^{d),Sj]. We now allow the system to be composed of n physical 
qudits. Define 

e,, = Y.^^^\ (88) 

i 

where the superscripts identify qudits p and g, all other qudits are acted upon by identity 
operators, and we have written the sum explicitly for emphasis. We can now show that 
all of the generalized exchange Hamiltonians Cpq for this larger system also commute 
with the collective errors Sa- To see this, note that the collective errors acting on the n 
qudits are again defined as 



n 



S^ = ^X^;^\ for a = l,2,...,rf2_i. (89) 



The only nontrivial contributions appearing in the calculation of [cp^, 5*^] have the form 

[A„A,](^)A;'^) + AfMA.,A,](^) (90) 

with the identity lying at all other positions. Since these terms may be rewritten as 

2^^-.(A?^Al^) + A?)Ai^)), (91) 

or 

^^k^i>^^>^ - A^ a!^)) = 0, (92) 

we see that the exchange Hamiltonians are also compatible with a DFS/NS that is 
supported by n physical qudits. 
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This generalization to n qudits is also applicable to the Hamiltonian F{d). Recall 
the commutator of F{d) with a three particle collective error. In this case we have 

[F{d), Si] = [f,jk\ ®Xj®Xk,\®t®l+l(^Xi(^l + l®l®Xi] 

= fijki[Xi, Xi] (g) Xj (g)Xk + Xi® [Xj, Xi] (S)Xk + Xi® Xj (g) [Xk, Xi]) 

= '^ifijkifilmXm ® Xj ® Xk + fjlmXi ® Xm CS) Afc + fklmXi ® Xj ® Am)- 

(93) 
This can reduced using Eq. IA.9[ After some relabeling we obtain 

[F{d), Si\ = 2i{djikdkmi — djikdkmi + djikdkmi 

djikdkmi + djmkdkli djmkdkli) Xi Aj Am 

= 0. (94) 

Now suppose that we have a system of n qudits with the interaction coupling particles 
p, q, and r, then 

[F{d), Si] = [f^jkX'fhfht^ , A, ® ... ® 1 + ... + 1 ® ... ® A,] . (95) 

Again, all of the terms in this expansion that do not contain ordinary products of two 
A's will cancel trivially, those which remain can be expressed as 

/...([A., A,](^) Xfxt^ + X? [A,, A,]('^) Xt^ + X^hf [Xk, Xi]'^\ (96) 

which becomes 

[F{d), Si] = 2i{djikdkmi — djikdkmi + djikdkmi 

— djikdkmi + djmkdkli — djmkdkli) Aj A^- A„ 

= 0. (97) 

The D{d) Hamiltonian can also be shown in a similar way to commute with the 
collective errors acting on a system of n qudits. Therefore, these DFS/NS compatible 
operators can be used to manipulate encoded quDits of arbitrary dimension as well. 

7. Conclusions 

We have presented a set of Hamiltonians that are compatible with an n-qudit DFS/NS 
[d and n arbitrary) encoding which protects against collective noise, i.e., any noise 
which affects the qudits in the same way. We have also provided a set of unitary 
transformations which can be used to perform any rotation on a logical qubit that 
is represented by the two degenerate representations in the product space of three 
qudits. The generalized exchange operations are analogs of the Heisenberg exchange 
transformation which is known to be universal for various systems of qubits and we 
have shown that they, as well as the logical gating operations, can be exponentiated 
analytically. In this context, we note that the three-qudit DFS/NS is similar in structure 
to the three qubit DFS/NS. 

Our analysis has focused on Hamiltonians which commute with elements of the 
stabilizer and the unitary transformations arising from the exponentiation of the 
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generalized exchange operation for qudits. For three qutrits we have found the complete 
set of Hamiltonians satisfying this commutation condition. However, these results could 
be extended by using the less restrictive requirement provided in Eq. ( l22ll . 

In order for these Hamiltonians to enable universal quantum computation they 
must also be able to generate entanglement between two encoded qubits. The CNOT 
gate acting on a subspace of the three-qubit DF-subsystem was provided in Ref. [38] 
using a circuit of 19 exchange interactions. Since the states considered there also appear 
in the three-qutrit DFS/NS as 1^2 and 1^2 we find that the exchange interaction alone 
can implement universal quantum computing over the three-qutrit DFS/NS as well. 
Furthermore, since the tableau for all thiee-qudit DFSs have an identical structure, 
those same states must also appear in the expansions of the logical states for d > 3. 
Therefore, the Hamiltonians derived here are sufficient to perform universal quantum 
computation using the three qudit DFS/NS. 

Finally we note that our algorithm for the determination of compatible 
Hamiltonians is also applicable to DFS/NSs which are not collective. In addition, the 
generalized exchange interaction is compatible with any collective DFS/NS and can 
be analytically exponentiated to produce the corresponding unitary transformations, 
including SWAP gates, for qudits encoding quDits. 
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Appendix A. The Algebra of SU{d) 

We have chosen the following convention for the normalization of the algebra of 
Hermitian matrices which represent generators of SU{d). 

tr(AiA,) = 2Sij. (A.l) 

The commutation and anti-commutation relations of the matrices representing the 

basis for the Lie algebra can be summarized by the following equation: 

2 

AjAj = -6ij + ifijk^k + dijkXk, (A-2) 

where here, and throughout this appendix, a sum over repeated indices is understood. 
As with any Lie algebra we have the Jacobi identity: 

filmfjkl + fjlmfkil + fklmfijl = 0. (A. 3) 

There is also a Jacobi-like identity, 

filmdjkl + fjlmdkil + fklmdijl = 0, (A. 4) 
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which was given by Macfarlane, et al. 

The following identities, also provided in [67], are useful in the derivation of the 
commutation relation given by Eq. (HSj), 

diik = 0, (A. 5) 

dijkfijk = 0, (A. 6) 

fijkfijk = d6ii, (A. 7) 

d^ -4 
d.,kdi,k = -^-Su, (A.8) 

and 

2 

lijmlklm 'J\^ik^jl ^il^^jk) ~r \0,ikmO'jlin Oijkm^^ilm) V-^^'^J 

and finally 

Jpiqjgjrjrkp I O^ / •'*i'^' (^A.iUJ 

(^piqjqjrjrkp \ o^ ; ij^i \^'^^) 

U'pigO'qjrJrkp I T^3 I Jijk-i {■^■^■^) 

/(P - 12\ 
O'piqO'qjrO'rkp I J73 I ^ijk- (^A.ioJ 

The proofs of these are fairly straight-forward and are omitted. 

Appendix B. Analytic Expressions for Unitary Transformations 

It was stated in 15.41 that the sum of all terms having the form Aj ® Aj can be shown to 
commute with the sum of all terms Aj ® Aj, where the Aj (Aj) represent the diagonal 
(off-diagonal) basis elements of SU{d). This appendix provides proof of this statement 
and gives an alternative description of the time evolution of the sum of off-diagonals. 

Appendix B.l. Sum of Diagonals and Sum of Off-diagonals Commute 
The objective of this appendix is to show that 



7 ,/^mo^ 7 ,/^jjc 



0, (B.l) 



where all Aj are diagonal, all Aj are off-diagonal, and ^lijk = \i ® Xj Cg) A^. It will 
be convenient to let i & I, where / represents a subset of all numbers in the set 
1, 2, ..., c?^ — 1 which correspond to matrices Aj which are diagonal and similarly j G J, 
with J corresponding to indices with Aj being off-diagonal matrices. 
First we define 

Cod = ^[Xi<» Aj, Aj (g) Aj] , (B.2) 
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and then proceed to show that C, 



od 



0. Here again the convention is that repeated 

we may rewrite 



indices are to be summed. Using the identities provided in Appendix A 
this last equation as 

Cod = — '^fijkfiji^k ® A/ + 2ifijkdijiXk 
+ "^fijkfijiM ® Afc — 2ifijkdijiXi 
= 2i{fijkdiji — fijidijk)Xk ® A/, 



)A, 
Afe 



(B.3) 



where the first equahty follows from Eq. flA.2p and the second from relabeling. Now, 
since fimk = for both indices i,m E I we may establish the following equahty 

/ J Jijk^^ijl / ^ Jirak^irali \^'^) 



«J 



i,m 



where m = l,2,...,(i^ — 1. Since the index i now serves to distinguish the diagonal 
elements, the Jacobi-like identity, Eq. ( 1A.4I) . can now be expressed as 

^ limk^^iml ~r Jlmk^iim U. yD.Oj 

This implies that if 

J2dni = 0, (B.6) 



then Cod = and the result follows. (Note that this is certainly true when i takes all 
values from 1 to d"^ — 1. Here, the i are in / so it may not be obvious.) To show this 
equation is satisfied for all d, we first state that it is true by direct computation for the 
matrix representation of the Lie algebra of SU(3) using the Gell-Mann basis and also 
for SU(4) using the Gell-Mann-like basis for SU(4). (It is trivially true for SU(1) and 
SU(2) since all dij^ = 0.) The only nonzero dm for z G / for SU(3) are 



d338 = l/v^5 and d; 



888 



-1/V3, 



and for SU(4) they are 

4,3,8 = I/VS, 4,8,8 = -I/VS, 4,3,15 = l/Ve, 

4,8,15 = i/Vg, 45,15,15 = -2/Vq. (B.7) 

Clearly, in both cases, J2i ^m = 0- We now proceed by induction, showing Eq. (IB. 61) to 
be true for d while assuming it is valid for d — 1. We also note that this assumption and 
proof is motivated quite well by the 4jfc for SU(3) and SU(4). 

By convention, we take a Gell-Mann basis for SU(d) which has its diagonal elements 
of the form 



Ad2. 



d{d 



1 ^ ° 


••• 





1 


••• 








1 ■■■ 






\ 



V 



(B.8) 



-{d-l) J 
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where there are d — 1 ones on the diagonal so that the matrix is traceless and normahzed 
such that Tr(Arf2_i)^=2. All diagonal matrices have this form and zeros are appended 
for higher dimensions. We now note that for all diagonal matrices Eq. ( JA.2P implies that 
the result of anti-commutation can only produce diagonal matrices. In other words, the 
Afc in Eq. (1A.2I) are diagonal whenever the left-hand side of that equation represents the 
product of diagonal matrices. Now, noting that 



{Aj, Xd'2-i} — 2a 



d{d-l) ^'' 



for z 7^ d^ — 1, we find that diia^-i = ^j2/d{d — 1) for all diagonal matrices Aj which 
are confined to the upper {d—1) x [d — 1) block. There are d — 2 such matrices since 
there are a total oi d — 1 diagonal matrices. The only thing left to find is the following 
dd?-i,d?-i,i- Let us calculate directly using Eq. (JB.SP and Eq. (]A.2p . 



{Arf2_l, Arf2_l} 



/I 
1 



d{d-l) 



^- 



V 

/ a 
a 
0a 

V ■■■ 








\ 



(rf-i)V 




h J 



where a = 4(2 — d)/d{d — 1) and b = A{d — 2)/d. This implies that the only nonzero 
dcp-i^cp-i^i is 

dcp-i4'2-i42_i = — (d — 2)^ 



d{d-iy 

and thus ^^ duj = for SU(d) assuming it is true for SU(d-l). D 



Appendix B.2. Exponential of the Sum of Off-diagonals 

Here we show how to calculate the exponential of ^ . Aj ® Xj, we will first show that 
the square oi K = ^Yljej^j ® ^j is diagonal. (It is not proportional to the identity; 
but it is diagonal with the same entry for each non-zero element on the diagonal.) 
To show this, we write the off-diagonal elements in one of the following forms, 

\k){l\ + \l){k\ , or t\k){l\-i\l){k\. (B.9) 

Note that K = J2k<i Qk,i so that K becomes 



K 



iE 



+ 



+ 



k<l 



+ ^ 



®(i 
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lJ2^i\m\ + m\) ® i\m\ + \m) 



2 

k<l 



= Y,[\m^\m + \m^\km- (b.io) 

k<l 

To square this, we calculate 

p<q,k<l 

X [|A;)(/| ® |0(A;| + |/)(A;| ® |A;)(/|] 

p<q,k<l 

+ Sqi6pk \q){l\ ® \p){k\] + 6pi6qk \q){k\ ® bX/|] 

However, note that the first term requires p < q = k < I and p = I. But p = I is 
impossible since p is strictly less than /. Therefore this term must vanish. Similarly for 
the last term. This gives 

K' = Y}\m\®\m\ + \m\®\m\)]. (b.h) 

k<l 

which is clearly diagonal. D 

In fact, the form is clear. This matrix has ones everywhere on the diagonal 
except at certain places. These places are when the indices coincide, i.e., when 
k = I, since we are restricted to have k strictly less than /. For example, using the 
Gell-Mann basis for SU(3), a straight-forward calculation using this formula produces 
ir2=diag{0,l,l,l,0,l,l,l,0}. 

It is now readily shown that 

K'^ = K, which implies K^ = K^, etc. (B.12) 

This allows us to analytically calculate the exponential of K and thus Cj which differs 
from K by the addition of diagonal terms (see Eq. fl83|) ) and tensor products of the 
identity, 

exp[-itK] = (1 - K^) + K'^cos{t) - iKsm{t), (B.13) 

and thus we obtain SWAP by a slightly different method. 

Appendix B. 3. Exponential of the Logical X Operator 

The objective is to calculate exp[— iaX] and exp[— iaZ] analytically. There are several 
ways that we could do this. Here we choose what we believe is the most straight-forward. 
Let us recall the definitions 

^=y^(ei-e2), (B.14) 

and 

Z=\{ei + e2-2e^), (B.15) 

o 
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where 

ei = ^Ao® Ai® Ai, (B.16) 

i 

62 = ^Ai® Ao® Ai, (B.17) 



and 



63 = ^Ai® A, ® Ao. (B.18) 



The diinensionahty of 1 should be clear from context and we will use 1 both for a d-state 
system, 1^ and also the identity of a composite system 1^ ® 1^ ® ■ ■ ■ ® 1^. 



The following identities may be shown, using the identities in Appendix A 



e^ = ±{d'-l)l-^ei. (B.19) 



Products of two have a cyclic property: 

2 
d 



2 
6162 = -^Aj®Ai®]l + i^ fijkXj ® Ai ® Afc 

ijk 



and 



+ ^ dijkXj ® Ai ® Afc 

ijk 

= ^63 -tF + D, (B.20) 

eie3 = -e2 + iF + D, (B.21) 

d 

6263 = 361 -iF + D. (B.22) 

d 

Finally, we note that 

f A\ /d^ -A\ 12 

where i = 1,2,3 and i ^ j ^ k ^ i. 

At this point the calculation proceeds in a straight-forward albeit tedious manner. 
One simply computes X^ = [(-^l/12)(ei — 62)]^ and Z'^ = [(1/6) (ci + 62 — 263)]^ using 
the identities in this appendix as well as those in Appendix A[ After showing Z^ = Z, 



we know that Z'^ = Z^ (and similarly for X) so the series may be summed to obtain 
the desired analytic expressions for the associated unitary transformations. 
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